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ABSTRACT 
The general voroperties of heat exchangers are examined in 

this study. Mathematical models are developed for the mechanism 
of heat transfer that take place during 

1. Evaporation in two-phase flow 

2. Condensation 

Sib | eleyeheatigies (ora (ekeye pial gifs: 
From these models, computer programs nave been developed for the 
preliminary designing of the following heat exchanger tynes: 


1. Counter-flow cooler with no change of state 


2. Counter-flow heater with no change of state 
3. Counter-flow condensor 

in Counter-flow evaboravor with two-phase flow 
5. Cross-flow heater 

6. Cross-flow cooler 


These programs are designed to be used in the examination 
of the size requirments for heat exchangers used in the opera- 
tion of thermodynamic cycles. The cycles of interest in this 
Study are ones proposed for the Droduction of useful vower from 


sources not conventionally used at Bresent. 
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ts INTRODUCTION 


In recent years, many proposals have been made for methods 
to produce power from sources not conventionally used. With the 
recent increase in the cost of fossil fuel energy sources, it 
can be expected that many such proposals will be forthcoming in 
the near future. 

One such proposal came from Professor Clarence Zener of 
Carnegie-Mellon University(1). This proposed evcle would pro- 
duce power utilizing the temperature gradient in the tropical 
areas of the oceans. To do this, a working fluid with a high 
vapour pressure is heated by the water at the surface of the 
ocean. In the tropics, this is around 25°C, The working fluid, 
as a vapor, then drives a turbine to produce power. It is then 
condensed back into a liquid by the deep ocean water at a temper- 
ature of 5°C. 

Hilbert and James Anderson(2) in 1966 proposed a cycle simi- 
lar to that of Professor Zener. It relies on the heating of the 
ocean surface water by the sun as the means for powering the sys- 
tem. In 1966, the cost of such a power source was deemed unrealis- 
tic, therefore the project was discontinued. In the light of 
the current shortage of energy, this provosal might warrent re- 
examination. 

Another such proposal came from a husband and wife team. 
Aden and Marjorie Meinel, optical physicists, proposed a differ- 


ent way of using the solar energy to produce unsable power(3). 





Their proposed cycle would use solar energy collectors to heat 
a liquid metal. This liquid metal would heat the working fluid 
which drives a turbine to produce power. 

The use of geothermal energy, heat from the earth's ine 
terior, as a power source has also been suggested(4). This 
type of heat source would be used in a manner similar to the 
warm surface water in the Zener proposal. 

All of these proposed cycles require similar RIOR 
The heat source is used to heat or evaporate the working fluid. 
For this, a heat exchanger is needed. The working fluid then 
drives a turbine and is cooled or condensed in another heat ex- 
changer. A pump or compressor is used to raise the pressure of 
the fluid and the cycle begins again. 

The feasibility of these and other such proposals must be 
examined. One phase of this examination should be to obtain an 
estimation of the size of the equipment necessary for the opera- 
tion of the cycle. A major item, both in the performance of the 
cycle and the size of the equipment, is the heat exchangers that 
are required. 

It is proposed that this study be an examination of the re- 
quired heat exchangers. This will be carried out by the develop- 
ment of a computer program for preliminary design purposes. 
General properties of heat exchangers will be examined. Mathe- 
matical models will be developed for evaporation, condensation, 
heating and cooling heat exchanger design. 


The heat exchangers needed for the cycle proposed by 


ao: 





Professor Zener will be designed in this study. The author 
feels that this will be useful for two major reasons. First, 
it will serve as an example of how this study and the computer 
programs developed can be used. A careful examination of the 
zener proposal can also be useful in determining if it is a 


feasible method of producing power. 





If. HEAT EXCHANGER PROPERTIES 


ae 1 Generel Properties 


2a ae Introduction 


A heat exchanger is a device used to transfer heat bet- 
ween two fluids that are sevarated by a wall. It is the form 
of the separating wall and the flow pattern of the fluids that 
Gistinguish between types of exchangers. This study will examine 
shell and tube heat exchangers and tubular cross flow exchangers. 
Fins are attached to the walls of some heat exchangers to increase 
the heat transfer. The use and effectiveness of fins will also 
be examined. 

Cost evaluation is an important problem in the heat ex-_ 
Cnenger design. A minimum cost, including intitial:=cost and the 
@est of pumping the fluids, is generally desired. The alterna- 
tives in cost are examined in this study. 

2.1.2 Shell and Tube Exchangers 

A shell and tube heat exchanger is made up of round tubes 
mounted in a cylindrical shell. The axis of the tubes are para- 
llel to the shell. Shell and tube exchangers are classified by 
the direction of flow of the fluids and by the method of reducing 
thermal stress between the tubes and the shell used (5 (page 18-34)). 
Selection of the type of exchanger to be used depends on a variety 
of properties including cost, pvressures and temperatures. 

meds tube-sHeet ExchanpereiSuune simplest to fabricate ana 


nas the least cost. In this type of heat exchanger the shell and 
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and tubes are rigidly held. Because of this method of connec~ 
tion the thermal stresses are the highest of any type of shell 
and tube exchanger. Another disadvantage of fixed-tube-sheet 
exchangers is the inability to clean the shell side of the ex- 
changer by mechanical means. 

To eliminate the high stresses of a fixed-tube-sheet ex~ 
changer, an exchanger with an internal floating head can be used. 
With this type of exchanger the tubes can be removed for clean- 
ing. The cost of an internal floating head exchanger is high 
and there is much more of a leaking problem than with a fixed 
exchanger. Fot the same job, the two types of exchangers would 
have the same size requirements and the design method is the 
same. 

The most commonly used shell and tube exchanger is the 
fixed-tube-sheet type(6). For this study it will be examined. 

A factor that has a large effect on the design of a heat 
exchanger is the choice of fluid to be passed through the tubes 
and that to be passed on the shell side. A number of factors 
go into the making of this decision: ( 5(Page 18-35) ). 


1. Cleanability. The shell is difficult to clean and 
therefore should have the cleanest fluid. 


2. . Pressure. It is much easier and cheaper to made 
tubes to withstand high pressure than the shells. 
High pressure fluids should nass in the tube. 


Sy ecempernacure, “High temperature fluids cause high 
stress. This stress can be more easily allowed 
for in the tubes. High temperature fluids should 
pass through the tubes. 


4, Quantity. The best design can be obtained when the 
fluid with the smallest flow rate is pvlaced in the shell. 


rae 





A counter heat flow exchanger has the flow of the two fluids 
in opposite directions. Figure 2.1-1. In most exchangers the 
flow is not purely counter flow. In the entrance and exit re- 
Ponca Ol the exchanger, some cross flow exists. If the exchanger 
isemore tnan twice as lone as it is wide, this cross flow region 
is small and can be disgarded. 

Parallel flow exchangers have the flow of both fluids in 
Hie same directions. See Figure 2.l-1. As with the counter 
Plow exchanger, the end regions have some cross flow which is 
mesregarded. Parallel flow exchangers require more heat trans-— 
fer surface area, for the same temperature rise, than counter 
flow. SeeFigure 2.l-2. Parallel flow also has the disadvan- 
Bace, compared to counterflow,,that the fluid being heated can 
only have a temperature rise equal to fifty percent of the hot 
meuv@ Lemperacvure at inlet. Counter flow exchangers can, in 
mae limiting case, have the fluid being heated reach the hot 
fluid inlet temperature(7). For these reasons parallel flow 
exchangers will not be further examined in this study. 

2.1.3 Cross Flow Exchangers 

In a Single pass cross flow heat exchanger, the two 
fluids move perpvendicular to each other. Figure 2.1-1. 
Cross flow generally requires more heat transfer area than coun- 
ter flow heat exchangers for the same temperature rise. Figure 
2.1-1. Cross flow units are used in areas of svecial application 
es the exchanger can be made compact. These exchangers are most 
widely used in gas to liquid exchangers with finned surfaces(8). 


(More discussion of fins appears in section 2.1.4). 
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Figure 2.1-1 
Flow Patterns of Heat Exchangers 
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Cold Fluid Temperature Rise as a Percentage of 
Hot Fluid Temperature Loss 


Figure 2,12 


Comparison of Required Heat Transfer Area 
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The cross flow heat exchangers examined in this study will 
be tubular exchangers. The heating or cooling fluid is passed 
through a tube and remains unmixed. The working fluid flows 


across the tubes and is allowed to mix. 


2.1.4 Finned Surfaces 

Fins are extensions of the tubes that increase the heat 
transfer surface. They are most effectively used when the heat 
transfer is between a liquid and a gas. A rough rule-of-thumb 
for the heat transfer surface area of the fin is that the sur- 
face area for the gas and liquid should be inversely propvortion- 
al to the ratio of respective heat transfer coefficient( 7(page 
182) ). This rule must be slightly modified because the effici- 
ency of a fin decreases with an increase in area. If, for ex- 
ample, an air-water heat exchanger was desired, the area of the 
air side should be 10 to 20 times as large as the water side. 
This is caused by the heat transfer coefficient of water being 
about 500 Btu/ hr ft °F to 20 Btu/ hr ft °F for air. To make 
up this need for different areas, fins would be used on the air 
side. 

The spacing of fins in a heat exchanger is generally de- 
termined by manufacturing limits. The height is limited by the 
reauired flow area. A spacing of from 6 to 15 fins per inch is 
the normal range(9). 

In the heat exchangers examined in this study, fins are 
not weed. Cost. 1cman amporvent factor. ~The fabrication costs 


of fins is much higher than for non-finned tubes. If the fin 
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meaenoy G€ither integrally or metallurgucally bonded to the tube, 
the fin efficiency may be adversely affected(1l). This requires 
the fins to be seddered or welded to the tubes. This is an ex- 
pensive fabrication process. Because fins decrease the flow 
area, they increase the required pumping power for the fluids, 
increasing the operational costs. 

The nature of the fluids used in this study also discourage 
the use of fins. In the proposed cycles, the heat menieion is 
from a fluid in one phase to a fluid of the same flow or to a 
two phase flow. Fin performance in both cases is less than ideal. 

If a detailed design of heat exchangers was going to be 
Berried out, the use of fins should be included. As this study 
is a preliminary design to determine feasibility, examination 


of fins is not needed. 


Mel,5. Cost of Construction 
The estimation of cost of a heat exchanger is a diffi- 

cult task and is often only a very rough estimate. Many rea- 
sons for this exist. Companies that make heat exchangers and 
components generally consider cost of both materials and fab- 
rication highly proprietry. They are, therefore, quite reluc- 
tent to release this data. What available cost data there is, 
is a number of years old by the time it is published. With 
mie current rate of inflation, data as little as six months old 
Zs inaccurate. 

Cost of a heat exchanger is sensitive to special require- 


ments. A small matter such as rigorous quality control can 
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easily double the cost of an exchanger(1l). The use of spe- 
cial materials, unusual size materials, or non-standard shapes 
Gan greatly change the cost of an exchanger. All of these face 
tors are difficult to include in a cost estimation. 

In this study a rough cost is examined. The values are 
from 1960(12), and are therefore auite inaccurate. They do serve 
to give a method of comparison between different shapes. If two 
exchangers are designed to do the same job, examining their com- 
paritive costs, the best can be selected, regardless of the ab- 
solute accuracy of the numbers themselves. 

Two cost factors are examined; the cost of the material 
maxing up the shell, and the tubes, and the cost of fabrication. 
The fabrication costs include header joint cost (welding), cost 
me Crilling of the header sheet, cost of cutting tubes to size 
emia pOSitioning of the tubes, The cost does not include inspec- 
tion, sales costs, and shop overhead(12). The material costs 
are determined on a base size heat exchanger. Corrections are 
made to adjust the base size to the particular heat exchanger 
designed. The two correction factors used are the Tube Material 
mrevowmana the lube Leneth Kaccor. The tube material correction 
takes into account the higher costs associated with the use of 
special materials. The correction for tube length decreases 


[mmenewleneth imnercases. This accounts for the economy of scale. 


Sl 





where: 


and 


and 


Cost = CW + CTB 
CW = CWV x N 
CWV 


3.0 if Dout is less than 4s in. 


CWV = 12/5 x (Dout - .5/12) + 3.0 if Dout is 
greater than 4 in. 


1, 
CTB = 15 -[(apea = 190)"]> TLF x AREA x TMF 
5 


a 
TLF 1.3 ~ (LENGTH - Oe 
10 
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NOMENCLATURE 


AT 
TLE 


TMF 


EuOScCripts: 
if 


oO 


Heat transfer surface area 

Cost of materials to build heat exchanger 
Cost of fabrication per tube 

Cost of fabrication of heat exchanger 
Outer diameter of tubes 

Heat transfer coefficient 

Thermal conductivity of material 
Length 

Wetted parameter 

Rate of heat transfer 

Temperature 

Change in temperature 

Tube length factor 

Tube material factor 


Cross sectional area 


With fin 


WEthouc. fin 


SiG= 


2.2 EVAPORATOR PROPERTIES 


mee,k Introduction 

Boiling heat transfer is the mode of heat transfer that occurs 
when a liquid changes phase to avapor. There are two regimes 
of boiling, film boiling and nucleate boiling. The form that 
boiling will take in a given circumstance depends on the derree 
of superheat of the liquid. Figure 2.2-1 shows typical pool 
boiling data(1). Pool boiling occurs when the heating surface 
is submerged beneath the surface of the liquid. In this study 
the region of interest is where Tw - Tsat is between 5 and 15 
peprees, in the nucleate boiling region. | 

Nucleation isthe formation, growth and motion of bubbles. 
mor nucleation to occur in a liguid, that liquid must be super- 
heated, its temperature raised above the saturation temperature. 
The nuclei is formed on a foreign object in the liquid. This 
object is usually a cavity on the heating wall but it can be 
suspended foreign matter with a non-wetted surface (4, page 6)% 
The bubble grows by heat conduction to the liquid-vapor inter- 
mace. the size it will become is based on a balance between 
buoyant and surface tension forces. 

When boiling occurs in a confined region with force convect- 
jon present, two-phase flow results. It is this type of boiling 
heat transfer that is discussed in this study. A mathematical 


model for two-phase evanoration is developed. 


-O)- 





A€Td OUT SuTwWo) uotzerpey 
SUTTTIOg WITT eTqeas 


atts 
STqeysuy pue Zutltog 
SZe9TONN Tetqaeg y 


QIVFJAIIUT © Ssty 
seTqqng sut{Tog ejeaTony 





SUT LTeg e3eeToOnNN 


UOTIOISLUOH VIITy 


au=a3/n9a “(v/>) Jon 


1.0 10 100 1000 10000 


ol 


Tw-Tsat CF 


Figure 2,2=-1 


Pool Boiling Data 


21 





> Two-Phase Evaporation 

Two-phase flow is a pattern that contains both liquid and 
vapor. It occurs when a flow is subject to forced convection 
in a confined space, a tube, annuli, parallel plates. The ob- 
jective of the designer of a two-phase evaporator is to determine 
the heat transfer and the pressure drop characteristics of this 


flow. From this data the size of the required heat exchanger 


mean be found, 


The determination of the desired characteristics can be a 
complex task as these two properties are counled thermodynami- 
cally(i9¥, page 47). <A change in the heat transfer changes the 
ouality of the flow and the flow pnattern. This changes the 
pressure drop. At the same time a pressure change alters the 
quality of the flow and therefore effects the heat transfer 
characteristics. 

The flow regimes change as heat is addedQio, page 14). 

Figure 2.2-2. The flow enters as a subcooled liquid. Boiling 

will take place at the walls even though the liaquid as a whole 

is below the saturation temperature. This forms a bubbly flow. 

me moderate velocities, greater than 3 feet per second, the 

bubbles will be fairly evenly distributed. When the liquid reaches 
mhne saturation point, bulk boiling will occur. As the quality in- 
creases the flow becomes annular with a thin layer of liquid on 

mae walls and a vapor core. In this vapor core will be liquid 
droplets that will decrease in amount and size as the ouality of 


eee Ltlow increases. At some point the heated surface will become 
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fry, This 1s called the burnout point. From here on the heat 
transfer mechanism changes and must be dealt with in a separate 
manner(5, page 13-35). 

With the complex flows and the coupling of the heat transfer 
and pressure drop properties, many of the solution techniques 
depend on a combination of analytical and experimental studies. 

A basis for any solution for a two-nvhase problem is the assump- 


tion that the flow is fully developed(i7, page 48). 


mec. 3 Mathematical Model 

Two models can be used for two-phase flow(q, page 24). The 
"homogenous! model considers the two phase flow as a single phase 
With mean fluid vroperties. This model is also known as the 
'friction factor! model. The'sSenarated flow' model considers 
the phases to be segregated into two streams, one liquid and one 
of vapor. If the velocity of each of the streams is considered 
to be equal, the ‘separated flow' model becomes the 'homogenous' 
model. As one would expect, the ‘homogenous! model is valid for 
the flow of low quality with the separated flow more applicable 
mor the higher quality annular flow. For the mass Flow expected 
in this study, G less than 5001b,/ft*sec, the separated flow 
model gives better agreement with experimental data(Qajo, page 47). 
For this mathematical model the ‘separated flow' model will be 
used. 

The following assumptions are made for the 'separated flow' 
model: (1) 


(overleaf) 


= 





ee 


1. The vapor velocity and liquid velocity are constant 
but are not necessarily equal. 


2. There is thermodynamic equilibrium between the 
phases. 


3. Empirical correlations can be used to relate the two 
phase friction multiplier #*% and the void fraction a 
to the independent variables of the flow. 

It is this third assumption that allows one to separate the 
calculation of the heat transfer and the pressure drop. The coup- 
. ling effect of the properties is taken into account by the empiri- 
cal correlations. 


To find the heat transfer in two-phase flow, the standard 


heat transfer equation q = A UAT is used. 
Where: qQ=w Neo 


Because there is a great deal of variation in the value of U 
over the range from a quality, x, of 0 for pure liquid to a quality 
of 1 for pure vapor, the heat exchanger is broken into segments 


of equal Axe 


To find U: 1.14 14 Dout - Din 


Uo hy Rw Kmat 


The last term is generally disregarded as Kmat 
Dout =- Din 
is much less than hy, and h,. 
Domi yecue s OCal Neatstransier of Ehe heating fluid, 
the following relationships are used: (Q2) 
For Re < 2100 


hy, = 
We D° 


alae 





NOume 00 =< Re <- 16000 


s116 ke Pry es 
hp = eee (Ree ale5) 
D 


For Re > 10000 


OCR eres Che. 


D »@ uyee 


ny = 


These relationshivs are emvirical equations based on results of 
many experiments. Other relationships to find h can be used. 
(1, page 192). 

Hercalculate hy, the Local heat transfer foefficient for 
evaporation, with two-phase flow is a bit more difficult. Chen (23) 
postulated that there are two mechanisms which take part in the 
heat transfer process for boiling in two-phase flow. One is the 
macroconveétive mechanism that is associated with bubble nuclea- 
tion and growth as takes place in pool boiling. Chen further 
postulated that these two mechanisms are additive to their contri- 
bution to the total heat transfer. Others have suggested this 
idea. Rohsenow(@4) in 1952 and Bambil1(Q5) in 1962 both concluded 
that the heat transfer effects of boiling and convective flow 
when occuring at the same time should be additive. 

For the contribution of the macroconvective mechanism, ChenQ3) 
Meit that the heat transfer coefficient h,,, , should be calcu- 


dated in a manner similar to forced convection flow. 


(equation overleaf). 


D5 = 





a K pe -8 4 
hmac 7 0.023 5 Re *® Pr 


This equation is modified to show the effect of the two-phases. 


h = 0.023 Re, -® pr, « KL 
mac = 9-023 Re, ry D 
where: Gw D ‘ol _ x) 
Rez Ne 
Hy 
PY = Cpl HL 


k1 


This leaves the only unknown as F. F is the Reynolds factor and 


is a function of Xz. 


where: 
F = (Re/Re,) °° 


and a -~ x)? 5 1 
O u 
Xtt = ieee ods 
x PT Hy 


For experimental data a curve of F vs 1/X¢~ is plotted. (Figure 2.2-3) 





An equation can be empirically derived. 


he = Oe S ds 4 aT 
X¢t Xet 


The contribution of the microconvective mechanism is derived 
from the analysis of Foster and Zuber(26) on pool boiling. Foster 
and Zuber found: 


(equation overleaf) 


ele 





LO 


10 





LO; 


1/X 44 
Reynolds Factor 


=). 


a0 





Nu, = 0.0015 (Re,)-°? Pr, -* 


Nb Lb 
ky 





r, is the radius of the bubble and is derived from bubble buoyancy 


and surface tension forces, 





Where: AT k Co.dg)> 6 225 
| mp = ame ft 
AP g AP 


Forster and dZuber in their analysis disregarded the fact that the 
degree of superheat is not constant across the liauid film on the 
wall material. For convective flow this is important as the 
temperature gradient across this boundary layer is made much 
steaper by the flow rate. 

Chen(Q22) then rewrote the Foster and Zuber equation to take 


this temperature gradient into account. 


0.00122 k,°”? BS path 625 2 
mic ~ ate RD a AT #* ap 75 § 


h 


5 is the supression factor and is pnlotted as a function of 
me = Re, F 4, 

The total heat transfer coefficient h,, is the sum of Nymac 
and hmic - This correlation appears to give better comparison 
with experimental data than any other correlations(G, page 126). 

The heat transfer coefficient h,, cannot be calculated in the 


above method for the region that is after the burnout point. 
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Figure 2,2~4 
Supression Factor 
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For this study it is assumed that burnout will occur at a quality 
of .95. For the low amount of Superheat expected this appears to 
be a reasonable value. For the segiment of the heat exchanger 
from quality of .95 to quality of 1, the following should be 

| used for h, ;: 


Hn O23 Dey We, Pe 


hy = 
2 8 
D Ly 


Once the value of U is determined, the required heat trans-— 
fer area can be found for each segment. 


wih Ax 
Area = pena te ae 
U AT 


Where: 

a tyalt > 2sat 
By assuming a geometry of the heat exchanger, the number of tubes 
and the diameter of tubes, the length of each segment can be 


calculated, 
Az = Area/nmD n 


This completes the heat transfer calculations. 

To calculate the frictional pressure drop in two-phase flow, 
the Lockhart-Martinelli-Nelson model(2!, 27) will be used. Although 
this model was derived for adiabatic flow, it is valid for heated 
Pipes(S, page 14 - 9). The heat exchanger is broken into the 
Same segment as for the heat transfer calculations. The pressure 


drops are then summed for the total. 





The relation of two-phase friction pressure drop to single- 


mnase friction pressure drop is: 


. 


(on a 
Az TPR Az “5 


Where or, ° can be found empirically as is plotted as a function 


Cte Or @Xp eet ) 


haf 
AZjmpp (42 jy 


Where $1* can be found empirically 


Then 


$7 = ov 
Xtt 
as = oy? x 6 [ap] 
(AZ) ppp (4Z) 470 
Where: ee 
42] . °09 G, “ uv’? if Re>2500 
RZ LO g Py 1D ice 
And 
a2] _ 128 Gy vy if Re < 2500 
AZ LO g Py, D2 
p 
By substitution, (7) can be determined. AP, is calculated by 
ek 
multiplying (7) by the length of the segment Az. The sum 
AMA 


2 2— 





100 


01 ex LAO 10 100 


Kee 


Figure 2.,2=-5 
Two-Phase Frictional Multiplier 





meevuie AP is the friction pressure drop in the evaporator. 

If this pressure drop is too large to be handled by the 
cycle in which the evaporator is to be placed, the geometry of 
the evaporator can be changed. Increasing the flow area should 
decrease the length and the pressure drop. 

This concludes the necessary mathematical model for two- 


phase flow evaporators. 


= 





NOMENCLATURE 


er 


AP 
Az 


aP 


:Q 


Re 


aT 


Heat transfer surface area 
Specific heat at constant pressure 
Diameter 

Inter diameter of tubes 
Outer diameter of tubes 
Reynolds number factor 
Mass flow velocity 
Gravitational constant 
Heat transfer coefficient 
Latent heat of evaporation 
Thermal conductivity 
Number of tubes 

Nusselt number 

Pressure 

Prandtl number 


Pressure drop per length 


Change in pressure 
Rate of heat transfer 
Radius 

Reynolds number 
Supression factor 
remperature 


Change in temperature 
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NOMENCLATURE (cont. ) 


Subscripts 
b 

h 

1 

1m 

mac 

mat 

mic 

sat 


cPr 


Quality 

Martinelli parameter 

Change in quality 

Mass flow rate 

Length 

Latent heat of vaporization 
Viscosity 

Density 

surface tension 


Two-phase frictional multilplier 


Bubble 

Heating fluid 
Pquad 

Lograthnic mean 
Macroconvective 
Material of walls 
Microconvective 
Saturation 
Two-phase flow 
Vapor 


Working fluid 
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2.3 CONDENSATION PROPERTIES 


mes... Introduction 
Condensation is the removal of heat from a vapor to convert 
it into a liquid. This is generally done by using a surface with 


its temperature, .T,, , lower than the saturation temperature, T 


W sat 3 


of the vapor. The vapor will condense into a liquid on the colder 
surface. The exact nature of this condensation mechanism is not 
fully established. It is felt that the condensation begins in 
omatl cavities in the surface much as bubbles are formed in bolil- 
ineg(S, page 12-2). 

There are two forms of condensation = dropwise and film. 
Dropwise condensation occurs on a non-wetting surface. Droplets 
moerm and run down the surface, increasing in size by further 
condensation and coalescence. New droplets are formed to take 
their place. The heat transfer for dropwise condensation is 
from 5 to 50 times that for film condensation(31). It is diffi- 
cult to maintain dropwise condensation for more than a short 
period of time without treating the heat transfer surface. 

Film condensation occurs when the liquid being formed wets 
the surface and establishes a stable film. This type of conden- 
sation is the prevalent in commercial condensers. As the heat 
eransfer rate is lower with fllm condensation, condens rs are 
designed for this form of condensation. If the condensing should 
be dropwise for part of the operation of the condenser, it will 


moerorm better than design. In this study, film condensation 
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will be examined. 


Seer Film Condensation 

The first attempt to analyse film condensation was carried 
out by Nusselt in 1916. When the stable film is formed, vapor 
condenses on the liquid. The condensation takes place at the 
liquid-vapor interface because the heat is being transferred 
through the liquid. At this interface there is a continual 
interchange of molecules being condensed and molecules evapora- 
ting, with the net flow condensing. 

In his analysis Nusselt made the following assumptions: 
Q0, page 314). 

te tne Tlow of condensate in the film is laminar. 

eo. The fluid properties are constant. 

3. The condensate temperature distribution is linear. 

4, There are no changes of momentum through the film. 


5. The vapor is stationary and there is no shear force 
at the ligquid-vavor interface. 


6. Heat transfer is by conduction only. 
These assumptions are important in the later development of the 
mathematical model for condensation. The only modification to 
these assumptions is in the temrerature distribution. The actual 


temperature distribution is slightly curved(10, pare 239). 


3.3 Mathematical Model 
In the design of a condenser, the desire is to determine the 


Bequired heat transfer area. The rate of heat transfer, a, is 
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found by: 


Where: 


As was done in 


disregarded. 


and 


mere AT), is 


q =U AAT) 


e l_ + 1 + Dout - Din 
hw Ae Kmat 


qc|~ 


boiling heat transfer, Kat’ pout -~ Din can be 


q W Dep 
_ (Tsat -Te1) - (Tsat - Te2) 
AT im nn 
in ae - Tod 
Tsat ~ Te2 


the ‘logarithic' mean temperature. 


To find the required area, all that needs to be calculated 


is the overall 


hneatetransier coeiiicient, Us. The heat transfer 


coefficient for the cooling fluid, h,, is found by the following 


relations(22). 


if Re, < 2100 


Races OU CRE nme ET z) “s 


Pe LUO Shere <2 10000 


k 
he = 116 — OR Ps hale) oe 


(cont. overleaf) 
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sine) IRS > aliovaye 


0.023 ke Pro’ Ga” 
la D2 


YS a 
Hower iVvestinec cOclticientsot "heac transier for condensayion, 
hy, one begins with the Nusselt analysis. For condensing on 


Severticle plate, 





Y 3 oad 
h, = 0.943 EP (P, - Py) ky" hee 
[yy AT 


Where hf, = hp, + .68 Cp AT. 


g 
The deriviation of this equation is examined in detail in References 
me, 3hand/c. 

iiecscoeciricilent must be modified for condensation on a hori- 
montcal tube. The relation above is valid with the modification 


Maaco £ Sin is used to replace g. The average value of h for the 


tube is found by integrating Dy from 0. to 





Doutyy aval 


For a number of tubes in a vertical row, all of the conden- 
sate dropping from any tube is assumed to fall on the next lower 
pube. This will increase the film thickness of lower tubes and 
mBequire a modification in h, . 

<cieCe 
~ P,) oT Dep 


ny Dout Vy AT 


Hep, (oe, 





h, = 0728 


SN i= 





Ny is the number of tubes in a vertical row. Observations by 
Grant and Oswent@G4) show that the condensate seldom falls from 
the upper tube to lower tubes as a continuous sheet. This would 
Bend to lessen the importance of having the np term. 

Myo Other effects can change the heat transfer coefficient. 
mmese are the presence of noncondensable gas and changes in the 
surface geometry. 

The presence of even a small quantity of non-condensable 
gas in a condensing vapor has a profound effect on the heat 
transfer coefficient. The non-condensable gas tends to form at 
the liquid-vavor interface. This increases the partical pressure 
of the gas at the interface and increases the resistance to heat 
transfer. The details of the methods to calculate this effect 
mee piven in Reference3}). In this study, the effects will be 
m@asregarded. in most commercial condensers, ejectors are used 
to remove these gases and alleviate the problem. — 

Changes in the surface geometry of the tubes can greatly 
mncrease the heat transfer coefficient with film condensation. 
One of the most promising of these geometry changes is fluting 
the tubes. This geometry was first examined by Gregorig(33). 
unm experiments Gregorig obtained an h for condensing of 
8000 Btu/hr. ft? °F, a factor of four above what is obtainable 
with round tubes. Because of the complex shape of these tubes, 
it is difficult to analytically evaluate the heat transfer co- 
mig icient. Instead, experimental data must be used. 


Upon finding ne and h,, the overall heat transfer coefficient 


sie 





can be found and the required heat transfer surface area 
determined. A heat exchanger geometry, tube diameter and number 
of tubes is ehosen by the designer. The condenser length can 


be found from this. 
1 = Area/tn Dout 


The frictional pressure drop of the cooling fluid can be 





"calculated. 
af Re. > 2500 
=. Gots BE Se 
SE Pee er 
ae, 
Re, Dees OP. 
Din Go 
Where Re, = 
7" 


ai Re. <a 2500 


5 eel MP aleaG' 
~ Re, Din 2g 0% 


If this pressure drop is too large for the desired apvlication, 
the condenser geometry, tube diameter and number of tubes, can 
be changed to decrease the pressure loss. 


This concludes the necessary mathematical model for con- 


mensing on the outside of tubes. Condensing on the inside of tubes 


is not often carried out as it involves two-phase flow and is diffi- 


cult to analyse. 


=o 





NOMENCLATURE 


Cp 
Din 


Dout 


nb 


rY 


AP . 


Re 


AT 


Subscripts 


Cc 


Heat transfer surface area 
Specific heat at constant pressure 
innersdilamever Of cubes 

Outer diameter of tubes 
Gravitational constant 

Mass flow rate w/Area 
NeaiGwematic te recoeCt eC LeMme 

Latent heat of condensation 
Thermal conductivity 

Length 

Number of tubes 

Number of tubes in vertical bank 
Prandtl number 

Change in pressure 

Rate of heat transfer 

Reynolds number 

Temperature 

Change in temperature 

Overall feat transter coeiiacient 
Flow rate’ 

Absolute viscosity 


Density 


Cool ne J 1usd 
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NOMENCLATURE (cont, ) 


1 Liquid 

im Logarithmic mean 
mat Material of the walls 

sat Saturation 

Vv Vapor 

W Working fluid. 


ee ee 





2.4 Cooling and Heating Properties 


eo4.1 ine roduction 

In this study, the evaluation of the pronerties of cooling 
a working, fluid and heating A OrEane fluid in a heat exchanger 
mre carried out at the same time. Thermodynamically the mechan- 
ism of heating and of cooling is the same. Any differences that 
do exist in the relationships will be discussed in the develop- 
ment of the mathematical model. 

As was noted in section 2.1, counter-flow heat exchangers 
and cross-flow heat exchangers will be examined in this study. 

The heat transfer mechanisms of both exchanger typves are similar 
but the approach to the designing of each is different. 

The job to be done by these heat exchangers will be to trans- 
fer a specific amount of heat from one fluid to another. If the 
exchanger is a heater, the working fluid will be the cooler fuses 
For a cooler, the working fluid must be the hotter fluid. The 
mass flow rate, w, of both fluids in each type of exchanger will 
be specified as will be the inlet and outlet temperatures. This 


will be true regardless of the flow vattern of the faa Sy 


4.2 Counter-Flow Heaters and Coolers 

The mathematical model for the Sizing of counter-flow 
heaters and coolers is based on the solving of similtaneous equa- 
mons. This is accomplished by setting groupvs of known constants 


equal to a large constant, Kl, K2, K3, etc. The heat transfer 


See 





area, number of tubes, and the exchanger length will be deter- 
mened in terms of these constants. This method is described in 
detail in Reference 7. 

Mor take this model valid fer both heating and cooling, the 
hotter fluid in the exchanger is designated by a subscripnt 1 and 
me cooler fluid by a subscript 2. In a heater the working fluid 
is colder and is designated by the subscript 2, while for a cooler 
the inverse is true. The mathematical relationships aperotnecwies 
maenticai. 

A limiting constraint on the design of a counter-flow ex- 
meanger is the pumping power required to overcome the friction 
pressure loss of the fluid flowing in both the tubes and the 
shell. In this design method, this constraint is set by estab- 


lishing an allowable pressure drop, AP. 


If Re, > 2500 


Oa0e2 u 0.2 CG 1.8 le 
AP, = eo Aas (35) 
gr Py Din 1.2 
e 1.8 9 
AP, = kK1 1 


Din 1.2 
rf Re. <e2500 
(128 Wa" \ (1,18 7, 
aL = | i | a 


ee od ae 


ic 
UG. (G5 PIO 


hha 


Bice 





an COME PLLUTe 
a2 1.8 
Ps «(Pees ve Gol*r, 


g P5 ; Deqoc 


1,8 
bs =e Mae 


Dea 1,2 
1m Re»5 <2500 
2 1.8 1.6 
py = (228 we Gy) ie Ce 
E05 Reo Deq?? Deq !? 
Where: Deq = (1.271 (s/Dout)* -~ 1) Dout (7, page 324) 


This is for square pitch exchangers with s the spacing between 


tube centers. 





Din (wy/ (1 em) 


Dea* ) 
4 


Where: Rey 


and Rey = Deq (wo/ (1 


Next the heat transfer coefficient, h, is calculated. (30, page 442) 


If Re, > 10000 




















eh | ‘ ie: «8 
h) = OO By See Nee) Gy = 1 
| uye Dirt? Dias 
ite 2 500. -< Rey < 10000 
, ry Pp 33 G48 14.8 
h, = 116 ki Pri™* G1 = x3 Gl 
uy" Re ° Din? Dir?’ 
If Re, < 2500 
1.86 ky Pri G4 G48 
h, = ear i = 43 21 
o.oo Din? Dine 


Ll i 
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If Rey > 10000 

















o& ,8 8 
ee eee = Kl 22 
Hee Deq? Deq? 
fi 2500 < Res < 10000 
33 : 
ho = ~116 ko Pro Go® = xy Go* 
8 13 2 Zz 
H, Re, Dea Deq 
Lag Res < 2500 
5 3 8 
ho = ie OG Ko Oe Go = Ky Go 
Ho? Res° Dea? Deq? 


Then: “> (both known ) 


: [>in] 


Ke |Go 


by ceftinition Gay AEBS a| 
Go Wo/AX9 


Where Ax is the cross sectional area of the flow. 


Boe Sa AN Walia ee | 

Axo = 1/4 Dea? n = 1/4 Deg Dout n 
By substitution 

Gy _ Wy Dout)* [xo K5|33 = x¢ 
Then Gy _ Wi Dea Pout . 

£ Wo Din? 


Sie 





From solving for AP j/APo and Gj/Go2 


iw oe ral Pea" 





ho K4 Go Din 
= K3 (y9 ma [Bows 1.2 j=2}" : K7 
K4 Kl Din od 


The overall heat transfer coefficient can be calculated 


by: 
il = i r ales + Dout =} Din 
U hy ho Kmat 


The last term can be disregarded as it is insignificant compared 
merchne other terms. This equation does not include the effect 


me U of scale deposits. For a vreliminary design it can be dis- 





regarded. 
1 1 al: 
—_ = 
Nig ie 
1 Le Daler 1 K6 Dea”? K8 Din*® 
a 
K3 Gy° K4 Gy Gy 
2 
Where ae 
K3°  K4 Ded? 
The heat transfer rate 
of a UN Me ~ 
Where heer Dan 
_ ee 2 
d= fo moan Ub A a ae Dit gl 1 ye. Gq Cp, AT, 


~4g- 





mere Aly, is the difference between the-inlet and outlet tempera- 
ture of the hotter fluid and ATy, is the logarithmetic mean 


temperature. 





Then ) 
L = Cree ua CEP Laie Gj Din 
4A T),, JU U 
= 1.2 1.2 
L = KO K8 Das CT 


then by substitution of 


i 
2 = —ee 
a1 Daenie Bote 
Into 
K 8 
AP, _ es Gy 
Dae 


Py t) 
G, = |\—————-|_ = kl 
Kl K8 K9 


By definition 


Q2 
4 
| 


= wWy/AX1 


Therefore Wy y 
ete aes 
Gy T Dy 


The value of Din and Dout are selected so L can be calculated. 


Koes bine Kao 


ee 
i 


and the heat transfer area can be found. / 


A = Dinn L 
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This concludes the necessary calculations for the heat exchanger 
design. This design method assumes that: 

(1) The specific heat of the fluid, Cp is constant. 

(2) The overall heat transfer coefficient is constant. 

If these factors are not constant, the design approach is 
Meer ied, the above procedures are carried out using average 
values of the fluid properties. A vanes of n and G) is determined. 
The exchanger is then divided into segments such that Cp and U 
Ban be assumed to be constant over the segiment. The design method 
is then repeated for each segiment with n and G] kept constant. 

The length is determined for each segiment. These lengths can then 


be added to get the overall length of the heat exchanger. 


2.4,3 Cross-Flow Heaters and Coolers 
In addition to the information about the job to be done 

by the heat exchanger noted in section 2.4.1, one must know the 
condition of the flow through the exchanger to size a cross-flow 
heater or cooler. (3, page 18 - 81). As stated in section 2.1.3, 
for this study it is assumed that the working fluid is flowing 
across the tubes and is mixed. This keens the working fiuid at 
a uniform temperature as it leaves the exchanger. 

In the development of the mathematical model, the equation 


me U AAT, is again used with a slight modification - 


q = UA FAT): 


(cont overlean. 


eee 





W : - - - 
here sn, - Gites ToD) (T, 4 We) 
me 
in “he. -= oul 
Th 7 To2 


AT mec is the logarithmetic mean temnerature for counterflow arrangfe- 
ment. This is not the integrated mean temperature A eS eel 
defined as: 


Fo® ATn/AT1 mo (30) 


end is a function of ee and Z. 


Where: Epetoe cea Mledeene  sentecyaVencss on Emesneal exchanger. 
355 Tao = Te1/Th1 = TA (30) 
and ) = 
5 i, eae 
ed = Mae 


F can be found on tables as a function of 4 ande,. F also 
Maries with the flow pattern, if one or both of the fluids are 
mixed. These tables can be found in References 5, 10 and 30. 

(iicmva Me sored can be found 270m ds=5W Cp, AT,. All of 


these values are known. To find the heat transfer surface area: 


A = q/U PAT, _. 


U can be calculated by the relations for h, and hj, examined 


1, 2 
mr the section 2.4.2 of this study. 


By selecting, an exchanger geometry, number of tubes and tube 


aio = 





Giameter, the length of the required exchanger can be determined. 
L = Area/nmDin n 


The frictional pressure loss is then calculated. The relation- 
ships are examined in detail in section 2.4.2 and need not be re- 
merived, If the pressure loss is too great, the geometry of the 
exchanger can be modified and the length recalculated. 

mois GeSiten methoa assumes Ghat U is constant throughout 
the length of the exchanger. If this is not a valid assumption, 
fne exchanger can be broken into segiments. The above design 
procedure is then used for each segiment with the sum of each 


segiment length the total length. 
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NOMENCLATURE 


Ax 
Cp 
Din 
Dout 


Deg 


= NOUN 


AP 
Pr 


Heat transfer surface area 

Cross sectional flow area 

SOeC Ie Meaveav Constance, pressure 
Inside diameter of tubes 

Outside diameter of tubes 
Equivalent diameter for shell side 
Mean temverature difference factor 
Gravitational constant 

Mass velocity 

Heat transfer coefficient 

Thermal COncuierlyiey 

Constant 

Length 

Number of tubes 

Frictional nressure drop 

Prandtl number 

Heavewrads Ler eta cc 

Reynolds number 

Tube spacing 

Temperature of cooler fluid 
memperavure of hotter fiuidc 

Log mean temperature difference 


Overall heat transfer coefficient 
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NOMENCLATURE (cont. ) 


subscripts 
1 
2 


mat 


Mass flow rate 

Temoerature coefficient 
Effectiveness of heat exchanger 
ViESCOS16y Of Liuid 


Density Or veld 


Hotter Tluid 
Cooler fluid 


Tube material 
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2.5 Material Selection 


ce). 1. Pnenrodue l 1on 
In the selection of material for the tubes and other 

heat transfer surfaces in the heat exchanger, it is desirable 
to have the heat transfer coefficient as high as possible. 
An increase in the heat transfer at the surface of the material 
produces a corresponding increase in U, the overall heat trans- 
fer coefficient. This will decrease the reauired heat transfer 
surface area as the area proportional to the inverse of U. 

This augmentation of the heat transfer can be accomvlished 
in many ways. The surface of the tubes can be roughened by a 
number of methods. The heat transfer surface can be left the 
meme ouc Che flow disturbed. This increases the heat transfer 
by displacement of the fluid or by setting up a vortex flow. 
mener methods of augmentation include vibration of the heat 
transfer surface, fluid vibration, and the use of electrostatic 
Fields. Of these methods of increasing the heat transfer, the 
best apvears to be the roughening of the surface material. The 
vibration and electrostatic methods require outside power to 


operate and are therefore less attractive. ( 5(page 10-2) ). 


Be Dec. Methods of Enhancement 

The first method for enhancing the heat transfer to be 
seriously considered, was increasing the surface roughness of 
the transfer material. The most successful for single-phase 
flow appears to be a method devised by Dipprey and Sabersky(14), 


and one developed by Kemeny and Cyphers(15). Dipprey and Sabersky 
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produced the roughness by electronically depositing copper on 
a rod coated with sand grains. The rod was dissolved leaving 
a tube of copper with close-packed sand grains on the inner 
surface. This method, for high Prandtl number flow, increases 
me nm, neat transfer coefficient, by a factor of two. Figure 
2.5-1. The Kemeny and Cypher method uses a helical groove in 
the material and a helical protrusion from the material. The 
grooved surface showed little improvement in the heat transfer. 
The helical protrusion did increase the h by up to a factor of 
two. Figure 2.5=2. As can be seen from Figures 2.5-1,2, the 
mearects of the enhanced surfaces decrease with the increase in 
the Reynolds number. These methods also show no advantage in 
two-phase flow. 

In two-phase flow, the heat transfer is increased by pro- 
moting nucleate boilling at as low a heat flux as possible. 

Young and Hummel(16) placed spots of teflon on the heated sur- 
face. This promoted nucleation as shown in Figure 2.5-3. This 
lowered the reauired level of superheat reouired to produce 
Boiling. With an increase of Reynolds number, the effect of the 
teflon spots decreases. 

The methods of increasing the heat transfer can also in- 
clude the use of a smooth surface with a higher k. All of these 
methods of increased heat transfer would be useful in the design 

of the heat exchangers for the thermal cycles..What is needed 
is a method of comparison of these increases heat transfer sur- 


faces to unenhanced surfaces. 
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Dippry and Sabersky Results 
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Figure 2.5-2 
Kenery and Cypher Results 
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Ae) Methods of Comparison 


A method of comparison has been developed by Professor 
W.M. Rohsenow and others working in the Heat Transfer Lab at 
M.1.17.(17). The data on the surfaces is presented based on the 


same heat transfer area and tube diameter of equivalent diameter. 
f Re* is plotted against j Re. Figure 2.5-3. 


Where: f Re*® = P/v 2¢p” De*a/A/V u? 


_ Ah Pr? Dea 
VY A/V Cpr 


and j 


The following are kept constant: 


8, P, Dea, A/V, uy Cp, Pr 


Therefore: 
f Re* ~ P/V 
and 
eevee AD or eu 
V V 
where P = WAP/p 


For the same mass flow rate and material properties: 


Ah 
Vv Vv 


For any flow arrangement there is one NTU vers e€ curve(18). 
Beere: _ ATrise in cold fluid 
Thot in - Teold in 
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hlaterial A 


f Re 


NTU 
Figure 2.5-4 


2h = 





The two materials can now be compared at two points. See 
figure 2.5-3. The mass flow rate, w, the Deq, A/V, Thot 4, and 
Teold in remain constant throughout this comvarison. 

At point A, tha values of w, V, and Re are the same as at 
point 0. This gives the same size and shape heat exchangers with 
different surface material. The amount of heat transfer and re- 
ouired pumping power will be different. These can be plotted as 
shown in Figure 2.5-5. 

Point B has the same mass flow rate, pumping power required 
and volume as point O. A heat exchanger operating at B will have 
the same volume as one at 0, but the shape changes and there is 
more heat transfer. This comparison is made in Figure 2.5-6. 

If the comparison of the materials is to be made while doing 
the same job, point C and point O are examined. At these vnoints 
the same heat transfer (NTU), mass flow rate and pumping nower 
occur. The difference is that less volume is needed with materi- 
al B (point C). Figure 2.5-7. 

Since A/V has been held constant, the area reauired at point 
C is also less. This can be an important comparison if the Hone 
exchanger needs to be as small as vossible. A cost comparison 
can easily be done. As the pumping power is the same, the only 


Bvariable cost is the material costs. 


Cost,/ Sauare foot x A 
Bos ‘one 





Cost ,/ Scuare foot ~* Ay 
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Figure 2.5-5 
Comparison of Exchanger at Point O and 


at Point A 
elie 
NLU, 
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Figure 2.5-6 


Comparison of Exchanger at Point O and 
allay icoyslie ar aa 


Se 





If Cr is greater than 1, A will be a cheaper material 
to use, If, however, Cr is less than 1, B is the cheaper 
material. 

A fourth comparison can be made at point D. Here the heat 
transfer, mass flow rate and the volume are kept constant with 
Geont O. The shapes differ and the pumping power is lower. 
Figure 2.5-8. If the cost of the reauired pumping power was 


too high with material A, material B could be used to lower it. 
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Figure 2.5-7 
Comparison of Exchanger at Point O and 
at Point C 
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Figure 2.5-8 


Comparison of Exchanger at Point O and 
at Point D 
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NOMENCLATURE 


or 


Deg 


NTU 


AP 


rr 


Subscripts: 
A 
B 


Heat transfer surface area 

Specific heat at constant pressure 
Cost) ratio 

Equivalent diameter 

Hine Clon nt aceor 

Gravitational acceleration 

tOecalmneavy seiconaremmcoch. 1Clent 
Thermal conductivity 

Number of exchanger heat transfer units 
Pumping power 

Change in pressure 

Prandtl number 

Rate of heat transfer 

Reynolds number 

Logarithmic mean temperature 

Volume of exchanger 

Mass flow rate 

Exchanger heat transfer effectiveness 
Mass density 


Absolute viscosity 


Material A 


Material B 
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meL.. RESULTS 


In this study, the objective use to develov a series of 
computer programs to carry out. preliminary design on heat ex- 
changers required for thermal cycles. For this develonvment, 
Meneral properties of heat exchangers were examined and mathe- 
matical models were developed for various modes of heat trans- 
fer. 

The computer programs were developed and tested. All of 
the programs performed as desired excepnt for the cost calcula- 
tions and the author feels that they should perform success- 
fully in studies of a broad range of inputs. The programs are 
explained and written in detail in the Avpnendix. 

Two of the programs were used in detail in an examination 
of the proposal made by Prof. Zener. See Avpendix G. The evapn- 
Oration and condensation heat exchanger Drograms carried out the 
desired preliminary desisn and showed the variety of designs 
that can be carried out with these vrograms. 

The Zener vorovosed cycle was examined to determine its 
feasibility. With the use of freon-21 as a working fluid, the 
evaporator must have a volume of two billion cubic feet for a 
100 lega-watt power plant. The condenser will have a minimum 
Baze of over three-hundred million cubic feet of volume. More 


detailed results are included in Avpendix G. 
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IV CONCLUSIONS 


The following conclusions can be drawn from the results of 


the study. 


iN. 


The computer can be a valuable aid in heat exchanger 
design work. 

The cycle vrovosed by Prof. Zener is not feasible at 

the present time. 

The cost analysis done in this study is very inaccurate. 
Although the programs can give a variety of geometrical 
arrangements there is no way to comvare these arrange- 


ments to determine tne best. 
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V RECOMMENDATIONS 


This study appears to the author to be weak in two areas. 
Tne evaluation of cost is very inaccurate. The problem of opti- 
mezation was not discussed in this studv. A further study of 
these two areas would be beneficial. As these two areas are 
inter-related, they could be examined together. 

More examination of the cost of heat exchangers needs to 
be earried out. The relationshivs for cost caluculation deri- 
ved in section 2.1 are based in 1960 data. This in in itself 
is a vroblem, but also the size of heat exchangers designed is 
way beyond the range of validity of these relationships. To 
extrapolate to 4,000 feet a curve that is valid up to !00 feet 
Seenrealistic. More ud to data values of cost should be de- 
Bermined. Also cost relationships that are valid for large 
heat exchangers should be derived. These relationships eourla 
easily be vlaced in the program. 

fiesprosrans as WriGlen CcCOntalm no ObtiImLZation operations. 
The programs reaulre a disigner to invut various peometric 
erransements and then compare the outnout to determine wich is 
eer fOr his application. With up to data cost relationships, 
my Optimization operation could be incorvorated into the pro- 
grams giving the designer the best heat exchanger within some 


mpecified limit. 
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Appendix A - COUNTER-FLOW CONDENSING HEAT EXCHANGER PROGRAM 


cr, 1 MSeceUuctlons fOr Obveration and input 

This program is written in FORTRAN IV and can be run as is 
on the Interdata computer of the Joint Civil and Mechanical 
Engineering Computer Facility at M.1.T. With modifications in 
the control, input and output cards, this vrogram can be run 
On any computer that is programed to comvile FORTRAN IV. 
| For this orogram, the inout data is broken into two grouns. 
mhe first contains four cards, three of which contain properties 
Bt the fluids and the fourth is the number, M, of input sets to 
follow. The second group has 2M cards. There are two cards to 
aset. The first card serves as an identification card and the 
second has the geometric arrangement of the heat exchanger; Din, 
ment, ©, N, KMAT, and NB. With this input peace mente M varia- 
tions of the geometry of the exchanger can be examined Ane the 
same fluid properties. 

The input data is punched on the cards as follows: 


Oe tees viddiepe ne unirer on elbe. hr., 
ECG el. mamGem er , 


(2) All values will be fixed point numbers with 
decimal points except the value of MM. 


(3) No commas are needed to separate the input 
values, 


(4) The input values must be in the proper col- 
Unitas dUGc meea moe be riehy or leit justified. 


(5) All fixed point numbers can have uv to 4 (four) 
CGEuGtnswalter the. decimal point. 


(6) The properties of the fluids, unless specified 
VOmDCeoueamOUarulCularewOunc., are average values 
evierytne length of the exchanper. 


ie 





CARD COLUMNS 


1--13 14—-—26 2{--39 40--52 53--65 66--78 
1. ex BEG 2 Tw | DELHW MKRW ROWL 
2 ROWV MUWL KWL KC PRC MUC 
1--12 13--24 25--36 37--48 49—-60 
. ROC HFG CEC ENG AL EG 
i: 1--4 
5 Plemr rit CAaGLenemame 9-1 Oreecacm feOMmermiCc Variation 
1--13 14--26 27--39 4QO--52 53--65 66--78 
6 DIN DOUT S N KMAT NB 
Be Same as Card 5 
8 Same as Card 6 


This repeats M times. 


The output will appear with the identification-at the top of 


mepage and the following printed below it. 


) MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE COOLING FLUID 
NUMBER Ob TUBES 
Penh DEAE ER OF suis 
OULD DIAMETER. OF 2UBES 
SPACTHG OF TUBES 
PENGTE OF TUBES 
DROP IN PRESSURE COOLING FLUID 
COMPARE Ly COSTS 


— 


For each identification there will be a separate page. 
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A.2 


AREA 
COST 


CPC 


CTB 


CW 
CWV 
DELHW 
DIN 


DLPC 


DOUT 
DT LM 
GC 
HC 


HFG 
HWA 


KC 


KWL 


KMAT 


LENGTH 
M 


LIST OF VARIABLES 


Heat transfer surface area in ft? 
Cost of heat exchanger in dollars 


opecific heat at constant pressure of cooling 
pital ala Jeng Wok ne 


Cost of the materials for the tubes and shell 
in dollars 


Nabmiecavul OnmCOsvonen GCOllars 
Fabrication cost per tube in dollars 
Change inenthalpy in Bin; 1b,, 

Inner diameter in feet 


Frictional pressure, Gropp in@eooling fluid in 
Tg V/ eee 


Outer diameter of tubes in feet 
Log mean temperature difference in OF 
Mass flow velocity in lb,/hr. ft.? 


NMeavee pans emucoct? telenucot cooling fluid win 
ISSO UV Apap ee 1D 


Latent heat of condensation in Btu/1lb, 


Average heat transfer coefficient of working 
fel a Ce eienie sey eines b> ee 


Inenmal sconcuecav@e ys omecooriat fluid in 
eV een ee 


racrmaleeondlecavicy OL worming fluid as a 
Ficus dn sei ar. oo. OF 


Peermal eonduerivacy (Of walle material in 
Bey litte eb, Ol: 


length of heat exchanger in feet 


Number of ~roups of inputs data 


le 








Mee LIST OF VARIABLES (cont. ) 


MFRC Mass flow rate of cooling fluid in IL oe gah ae 
MF RW Mass flow rate of working fluid in 1b,/hr. 
MUC PauscOchty  OlscoO hn ehlnece in, Lie iit. fit. 
MUWL VascOsmeyeOUmvOMmuime —limares ao .fquid in 
LS ARO 1G « 
N Number of tubes 
NB Number of tubes in a verticle row 
PCI Inlet pressure of cooling fluid in lb/ft.? 
PC2 Outlet pressure of cooling fluid in lb/ft? 
PRC Prandtl number of cooling fluid 
Q Rate of heat transfer in Btu/hr. 
REC Reynolds number of cooling fluid 
ROC Density sot scooltineuluideeian lb/ft. 
ROWL Density of working fluid as a liauid in TD Ges 
ROWV Density of working fluid as a vapor in lby/ft.°? 
S spacing of tubes | 
rC1 Inlet temverature of cooling fluid in OF 
rC2 Outlet temperature of cooling fluid in °F 
TW Temnerature of working fluid in OF 
TLF . tubo lonPenetaerLor in do) larc/it. 
TMF Lube Malverne taccox 
U Overall heat transfer coefficient in BTU/hr.ft? OF. 
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Appendix B -= COUNTER-FLOW EVAPORATION HEAT EXCHANGER PROGRAM 


=, 1 Inscructaons [or operaw:on and innut 

This program is designed to be run on the Interdata com- 
puter of the Joint Civil and Mechanical Engineering Computer 
Sacility at M.1.T. The program is written in FORTRAN IV and 
can be run on any computer that is vrogramed to comnile FOR-=- 
TRAN IV. If run on a comouter other than tPacmlnGerdata. some 
modification may be required to the control, input and outout 
Bards. 

The input data is broken into two sections. The first con- 
mains five cards, four contain fluid properties and the fifth 
me the number, M, of inout sets to follow. The second section 
Has 2M cards. The first card serves as an identification card 
and the second card has the geometric arrangement of the heat 
meechanger; DIN, DOUT, S, N, KMAT. With this inout arrangement, 
M Pie ferent variations of feometry can be examined using the 
same fluid pronerties. 

The input data is punched on cards with the following 
rormat: 


(SAID values will be in units of Ib), hr, 
WDleten ENG We tehelom CIM. 


Mejor values except will bertaxed point 
numbers and will have decimal voints. 
Hees) at LOAatIng “DOUNntL. number, 


(3) Wo commas are needed to separate the invut 
Values on a2 Single card. 


(1) The input values must be in the desirnated 


Clint Olt need nor De riphtG sor lert justi= 
fied. 
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(5) All fixed point numbers can have up 
to four digius santer the deeimal 
polinga 


(6) The fluid nroperties, unless snpeci- 
fied CO De saAr va particular woine, 
are average values. 


CARD COLUMN 
1--14 15--28 29—-—42 4 3—_—56 51--Te 
al DELHW selil ie TSAT ROWL 
2 ROWV MUWL MUVW CPWL KWL 
3 MPRW esol PHe PWi PWe 
4 HFG CPH MUH ROH KH 
1-—4 
5 i 
6 Identification name is columns 5----- 18 
1--14 
27 DIN DOUT S N KMAT 
8 Same as card 6 
9 Same as card 7 


This is repaeated MM times 


The output will appear with the identification name at the 
mop of the vage and the following vorinted below it. 


MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE HEATING FLUID 
NUMBER OF TUBES 

Pee DIAMETER OF TUBES 
OUTER DIAMETER OF TUBES 
SPACIWG OF TURES 


(cont. overleaf) 


—0= 





LENGTH OF TUBES 

DROP IN PRESSURE WORKING FLUID 
DROP IN PRESSURE HEATING FLUID 
COMPARITIVE COSTS 


Por each identification there will be a separate page. 


Se 





B.2 LIST OF VARIABLES 


AREA 
COST 


CPH 
CPWL 
f.B 


CW 
CWV 
DELHW 
DELPH 


DELPW (1) 


BELZ (1) 
DELX 

DEQ. 

DIN 

DOUT 
FXTT 


GH 
GV] 
HH 


HFG 


Heat transfer surface area in ft? 
Cost of heat exchanger in dollars 


SOeCCHLuUCweAabeae scomstane pressure of heating 
I IVE] Sigh Seu es S18 


Specific heat at constant pressure of working 
fluid as a liquid in Btu/lb,, °F 


Cost of the material for the tubes and shell 
in dollars 


Costrol “fabrication in dollars 
Fabrication cost per tube in dollars 
Cicnpe  iecuna pyar Vorine wlUid ta Btu/ ibe 


Change Anvpressurevos neating fluta In Ib/ft- 


Frictional wressure drop of working, fluid 
per segiment of exchanrer in 1b/ft? 


Length of each segiment in feet. 
Change in quality in each segiment 
Equivalent diameter in feet 

Inner diameter of tubes in feet 
Outer diameter of tubes in feet 
Reynolds number factor 


Mass flow velocity of heating fluid in 
lb/hr. ft? 


Mass flow velocity of working fluid in 
IS EV AGS ta te 


Meat transrer weoetiicient ot heating fluid in 
Beu/hr. £t* OR 


Latent heat of evavoration in Btu/1b,, 


209 = 





me2 LIST OF VARIABLES (cont. ) 


HMAC 


HW 


IV 
KH 


KWL 


M 

MF RH 
MFRW 
MUH 


MUWL 


MUVW 


N 
eid 
PH? 
PW 
PWe 
PRWL 
REWL 
ROH 


ROWL 


Heat transfer coefficient for macroconventive 
Such im Stuy newt. 


Heat transfer coefficient for working fluid 
chs) 2 Siejelene ois  SkeblAeban ibmm a 015 


Two phase frictional multiplier 


Thermal conductivity of heating fluid in 
Beuysnis tren 


Inermal COnguectivicy of Working fluid as 4 
ulslieithatiel alia Jeigul/ben ca4 210 


Number of groups of input data 

Mass flow rate of heating fluid in lbn/hr. 
Mass flow rate of working fluid in lb,/hr. 
Viscosity of heating fluid in lb,/hr. ft. 


MESCOsty tom Vorikiinmor luted as 4 11quid tn 
a 5/ Pee 


Viscosity of working fluid as a vanor in 
eyes a he 5 


Number of tubes in exchanprer 

Inlet pressure of heating fluid in 1b/ft*. 
Outlet pressure of heating fluid in lb/ft?. 
ileterrescsireros working fiuid in ib/ft-. 
Outlet pressure of working fluid in 1b/ft?. 
Peade lenumbper OUeworkKingstiluid as a liquid 
Reynolds number of working fluid as a liquid 
Density of heating fluid in 1lb,/ft?. 


Density of working fluid as a liquid in 1b,/ft?®. 
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Be 2 


ROWV 


REH 


TH 
TH2 
TLF 
TMF 


TSAT 


XTT 


LIST OF VARIASLES (cont). 


Density of working fluid as a vapor in 
ihe ae 


Reynolds number of heating fluid 

Spacing between tubes 

Inlet temperature of heating fluid in °F 

Outlet temrerature of heating fluid in° F 

Tube length factor 

Tube material factor 

Saturation temperature in Op 

Quality 

Martinelli parameter 

Overall heat transfer coefficient in Btu/hr. ft? CF 


Length of heat exchanger in ft. 
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Appendix C. COUNTER FLOW HEATER HEAT EXCHANGER PROGRAM 


Pi Instructions for Operation and Input 

This program is written in FORTRAN IV and can be run as is 
on the Interdata computer of the Joint Civil and Mechanical 
engineering Computer Facility at M.I.T. The program can be run 
on any computer facility that is designed with the ability to com- 
pile FORTRAN IV. If run on a computer other than the Interdata, 
some modifications may be required to the control, input and out- 
But cards. 

The input data is broken into two sections. The first sec- 
tion contains eight input data cards. The first seven of these 
cards containthe fluid properties and the eight is the number, 

BEC Cmertigonu sets in the second section. @ Ihe second section has 
Meecoras., The first card in each of these M sets serves as an 
identification card. The second card has the data on the geometric 
arrangement of the heat exchangers; DIN, DOUT, and S. With this 
input arrangement, M different geometrical arrangements can be ex- 
amined using the same fluid properties. 

ines input data is punched on cards as follows: 

C1) GILT SacI Y STS), SMpib Ides exe alia) (yon tide! Che eee obain e 
Regier Gu onic 4. 

(2) All values except “ will be fixed point 
numbers and will have a decimal point. Mis 
BweEweavine enorme numoer. 


(3) No commas are needed to separate the input 
dataron a single card. 


es 








(4) 


(>) 


The input values must be in the desig- 
nated column but need not be rirht or left 
jusvified. 


All fixed point numbers can have up to 4 
(four) digits after the decimal point. 


(€) The fluid pronerties, unless svecified 
to be at a narticular noint, are aver- 
efe Values, | Any Droperty that the Nase 
figure of its variable name is a number 
iSespeeii ied To. de -eat That -tennerature. 
FARD =: COLUMN 
1--12 13--24 25--36 37--48 49--60 
1 MEFRW nea ihe JEG/aL TWe 
2 PHI PH2 ey PWe MUH 
3 MUW Crh CPW ROH ROW 
4 KH KW DCPW1 PEE 2 DCPW3 
1--10 11--20 21--30 31--40 41--50 51--60 61--70 71--80 
5 WCPH1 DCPH2 DOPE TW3 aS DMUW1 DMUW2 DMUW3 
6 DMUH1 DMUH2 DMUH3 DKW1 DKW2 DKW3 DKH1 DKH2 
7 DKH3 DROW] DROW2 DEROWS DROH1 DROH2 DROH3 
1--4 
8 M 
5-------- 18 
9 Identification name 
1--12 13--24 25--36 
10 DLN DOUT S 
3! Same as 9 
12 Same as 10 
This repeats M times 
(cont. overleaf) 
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| 


The output will appear with the identification name at the 


top of a page and the following nrinted below it. 


MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE HEATING FLUID 
NUMBER OF TUBES 

INNER DIAMETER OF TUBES 

OUTER DIAMETER OF TUBES 
SPACING OF TUBES 

LENGTH OF TUBES 

DROP IN PRESSURE WORKING FLUID 
DROP IN PRESSURE HEATING FLUID 
COMPARITIVE COSTS 


Hor each identification there will be a separate page. 
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LIST OF VARIABLES 


AREA 
COST 


CPH 


CPW 


CTB 
CW 
CWV 
DIN 
DOUT 
DEQ 


DPH 


KW 


LENGTH 


Heat transter area in ft-. 
Cost of heat exchanger in dollars 


Specific heat at constant voressure of heating 
five) an Gy oe 


Ssvecific heat at constant pressure of working 
PALUELo aBolverepy/ alee, 940 


Cost of materials for shell and tubes in dollars 
BoD racarMOnEcOsue. aldol airs 

Fabrication costs per tube in dollars 

Inner diameter of tubes in ft. 

Outer diameter of tubes in ft. 

Equivalent diameter in ft. 


Frictional nressure drop in heating fluid in 
ly 2G aa 


frietional pressure drop of working fluid in 
ible Ate 3 


Log mean temperature difference in OF 


hROnper uy at temperature 1 
Property at temverature 2 
Proverty at temperature 3 


Mass fllow velocity of heating fluid in lb_/hr mie 
Masse. Lome Oct wmoumometimae | fla im 1b.) / ime oe 


thermalveonductavicy ei heating I luida in 
UY Aan ea ak 18 em no) 


iherinal conauctlvity of working Fluid in 
IS NGIDY Age comes Me mss) o 


Length of heat exchanger in ft. 


50 





LIST OF VARIABLES 


MF RH 
MFRW 
MUH 
MUW 
M 

N 
eal 
PH2 
PW) 
PWe 


PART 


THI 
THe 
TH3 
Ty 1 
TWe2 


TW 3 


(cont). 


acs yi lonmravesof heatania ilutd sine b./hr 
Mass flow rate of working fluid in 1Jb,/hr- 
ve ES eoctiy sOlenea tng iluid an ab 7h. ft. 
Was Costin ome votineu ml peaeen 1D n/ line ft. 
Number of input data groups 

Number of tubes 

inlet pressure of heating fluid in lb/ft. 
OuULletworessurer on eatin —lutd an Ib/it- . 
Inlet pressure of wforking fluid in lb/ft". 
Outlet pressure of Workine fluad in lb/Eft-. 
Temperature ratio 

Lrandti, numberref heating tiluid 

Prana. lLe number Ob WwOrtene on futd 

Reynolds number of heating fluid 

Reynolds number of working fluid 

Density of heating fluid in 1b,/ft?. 
Deris#ty son Vorking fluidetne lb/ft *. 

Tube spacing 

Inlet temnerature of heating fluid in °F 
Outlet temnerature of heating fluid in °F 
Intermediate tennerature of heating fluid in On 
Inlet temperature of working fluid in °F 
Owtlet temnerature of working fluid in OF 


Intermediate temperature of working, fluid in oF 


=op= 





LIST OF VARIABLES 


TW1 
TW 
LV 3 


TLI 


_ = x9 
fe —--XX12 


Coon a 


Inlet temperature of working fluid in F 
Outlet temperature of working fluid in F 
Intermediate temverature of working fluid in 
Tube length factor 

Tube material factor 

are grouned constants 

are grouped constants 


Above properties for each segment of exchanrer 
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Appendix D. COUNTER-FLOW COOLER HEAT EXCHANGER PROGRAM 


m, 1 Instructions for Oneration and Innuts 
ines program is written in FORTRAN IV and can be run as is 
on the Interdata computer of the Joint Civil and Mechanical ‘Engineering 
Computer Facility at M.J.T. The program can be run on any com- 
eecer facility that is designed to compile FORTRAN IV. If run 
Bn a computer other than the Interdata, some modifications may 
Mmmbeauarea to the control, imeut and output cards. 
iWeminouy deck @S broken Into two Sections. The first con= 
tains eight input cards. The first seven cards contain the fluid 
Broperties and the eighth is the number, M, of input sets in the 
mreOmansecti1on. Tne second section has 2M cards. The first card 
Bf each set serves as an identification card. The second card con- 
tains the data on the geometric arrangement of the heat exchanrer; 
fim, DOUT, and Ss. With this arrangement of the input data M differ- 
ent geometrical arrangement can be examined using the same fluid 
properties. 
The inout data is punched on cards as follows: 
GO wale val vece Vl bic ml meun i (scmicme lb ae Gi. 
feet, Etu and °F 
(2) All values except M will be fixed point 
numbers and will have a decimal voint. 
Veiswan floating point, number 


(3) No commas are needed to senarate the innut 
data on a single card 


(4) The input values must be in the desifnated 


Collum DULL GOenOG Mmeea torbe rwiphtvsor left 
Hust ited 


-100- 





(5) All fixed point numbers can have up to 
WMGhoUr) Cdieits after thesdecimad fount. 


(6) The fluid vroperties, unless svecified 
LO De thay a0 a Particular pein. are 
average values. Any pnronerty that the last 
PiCUre Ol 16S Versaule Mame Mora. numbereis 
Specified to be at that temperature. 

FARD COLUMN 
1--12 13--24 25--36 37--48 49--60 

l TCl TC2 Tw TWe PW 

2 Ev 2 eal Re2 MFRW MUW 

3 CPW MUC CPC KW KC 

4 ROW ROC 1a DCPW2 DCPW3 


M-10 11--20 21--30 31--40 41-+-50 512--60 61---70 71--380 


5 mercl DCPC2 HCE DS TW3 nes MONE DMUWe DMUW3 
6 DMUC1lL DMUC2 DMUC3 DKW1 DKWe DKW3 DKC1 DKCe 
7 DKC3 DROW] DROW2 DROW 3 DROC1 DROW2 DROW 3 

j---4 
; 1 

5-------- 18 
9 Identification name 


H=-12 13--24 25-—-36 
10 DIN DOUT 5 

i} Same as 9 

12 Same as 10 


nese last two cards are repeated M times 


(cont. overleaf) 
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The output will appear with the identification name at 


the top of a page and the following vrinted below it. 


MASS FLOW RATE WORKING FLUID 
MASS FLOW RATE COOLING FLUID 
NUMBER OF TUBES 

INNER DIAMETER OF TUBES 

OUTER DIAMETER OF TUBES 
SEC GeOnmUU BES 

LENGTH OF TUBES 

DROP PRESSURE SYCREING BLUID 
DROP IN PRESSURE COOLING FLUID 
COMPAKI Rive Cosa 


For each identification name there will be a separate page. 
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LIST OF VARIABLES 


AREA 


CPC 


CPW 


COST 
or 
CW 
CWV 
DIN 
POUT 
DEQ 


DPC 


DPW 


DTLM 


GW 


KC 


KW 


LENGTH 


Heat transfer area in ft?. 


ppeciiic heat at Gonstant oressure for cocling 
Qololicl siiqy engi see) 9 


Svecific heat at constant pressure of working 
fluid in Btu/lb,, ° 


Cost of heat exchanger in dollars 

Cost of materials for shell and tubes in dollars 
Fabrication costs in dollars 

Fabrication costs per tube in dollars 

Inner diameter of tubes init. 

Outer diameter of tubes in ft. 

Equivalent diameter in ft. 


ricle GlOnaimapreessure- arepeer COOlIne Tluad in 
By sat. 


Frictional pressure drop of working fluid in 
Moyne 


Loge mean temverature in On 


PrOpercy at temperature 1l 
Property at temnerature 2 
Property at temperature 3 


Macc lov sveloeity soleieOodime fluid in lby/nr. £53 
Mass flow velocity of working fluid in 1b,/hr. ft*. 


Thermal scCOndUuecuiva PyYaO1scOOoline fluid in 
Buu/ hie. wit ee oo 


Thermal conductivity of working fluid in 
Btu / Mires ecsuo hn 


bencth ol neat “exchanger in ft. 
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LIST OF VARIABLES 


M 
MFRC 
MF RW 
MUC 
MUW 
N 
ECL 


PC2 


TCL 
TC2 
aC 3 
TLF 


TMF 


(Gone ic 


Number of inputs data groups 

Vas Selous meavenon scoGlingaeluid aim be / hr. 
Masiseloverd ve so! worsens ete ras ine) /lan, 
Vols @re shine eal “@eodla aye shoutal aya Wek Merah kie 
KEISE Os! (yon auoricing (flurde in lb/ft. 
Number of tubes in exchanger 

Inlet pressure of cooling fluid in 1b/ft?. 
Outlet pressure of cooling fluid in lb/ft?. 
Inlet pressure on vorkinewuimid in lb/ft>, 
Outlet pressure of working fluid in lb/ft”. 
Prandtl number of cooling fluid 

Prandtl number of working fluid 
Temperature ratio 

Reynolds number of cooling fluid 

Reynolds number of working fluid 

Density of cooling fluid in 1lb,/ft?. 
Density Of Worcine tf louta inl b,/ ft. 
Spacing of tubes 

Inlet temverature of cooling fluid in °F 
Outlet temnerature of cooling fluid in OF 
Intermediate temperature of cooling fluid in he 
Tube lLlenztn facvoer 


Tube material factor 
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Hrot OF VARIABLES (cont). 


TW1 Inlet temperature of working fluid in °F 

TW2 Outlet temnerature of working fluid in °F 

TW3 Intermediate temperature of working fluid in oR 
Xl----- KO are grouned constants 

XX1----XX12 are frouped constants 

X------ (J) Above properties for each segment of exchanger 
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Appendix E. CROSS-FLOW HEATER HEAT EXCHANGER PROGRAM 


Bl imstruction tor Operattonvand Inputs 

Written in FORTRAN IV, this program is designed to be run 
on the Interdata computer of the Joint Civil and Mechanical 
Magineering Computer Facility at M.1.T. With modifications in 
meeecontrol, input and output cards, this program can be run on 
any computer that is designed to compile FORTRAN IV. 

The first three cards of the input data are used for proper- 
ties of the fluids. Ten cards are then needed to input the ten 
Be ten Maciek ohne nest card Contains 2 sancle number, 1. 

This is the number of geometric arrangements that will be examined. 
whe input cards that follow are in two card sets. The first of 
mhese contains an identification name, the second has the data on 
geometric arrangement, DIN, DOUT, S, N, and KMAT. With this in- 
put arrangement, " different geometric arrangements can be ex- 
amined with the same fluid properties. 

new input data wis punched tonmthe cerds as follows; 


(1) All values will be in the units of 


peices leet ep Cle wands a 


(2) All values excent M will be fixed voint 
numbers and will have a decimal point. 
Mis a tileating pointe number: 


(3) No commas are needed to senarate the 
ipwe CaGa Ona si ntlescarae, 


(4) The input values must be in the desig- 


nated columns but do not need to be 
riciueor eft Us tad eC. 
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(5) All fixed point numbers can have up to 
4 (four) digits after the decimal point. 


(6) The fluid properties are average values 
of these vroperties over the temperature 


range. 
CARD COLUMN 
1--12 13--24 25--36 37--48 49--60 61--70 
1 MFRW TWH TC AUsiet THC PH1 
2 PZ PW1 PWe CPW CPH MUW 
3 MUH KW KH ROH ROW 


ards 4 through 13 contain values of the matrix fF 


Theses values are punched 10 to a card with six columns per value 


1--4 
14 M 

5-------- 18 
15 Identification name 

1--12 13--24 25--36 37--48 49--60 
16 DIV Dour S N KMAT 
iW Same ds 15 
18 Same as 16 


These last two cards are repeated M times 


The output will appear with the identification name at the 


mop Of a page and the following printed below it: 


(cont. overleaf). 


Sle 





MASS FLOW RATE WORKING FLUID 
MASS PLOW RABE HEATING FLUID 
NUMBER OF TUBES 

IWNER DIAN ETER OF TUBES 

CURE RSP wena Or ee VUBhS 

Sree iG Oa huis mes 

CENG LEO SV USES 

DROP eee tro ine OmblniG: shia Weep 
DROP SLEReRESSURE QIBATING PEUID 
COMPARITIVE COSTS 


For each identification name there will be a separate page 


=a 





fm2 List 


OF VARIABLES 


AREA 
COST 


CPH 
CPW 


CTL 
CW 
CWV 
DIN 
DOUT 
WEP 
% 

EC 

1g 

GH 
GW 


HH 
HW 
KH 
KW 


LENGTH 
i! 
MF REI 


*DTLM 


Heat transfer area in ft”. 
Cost of exchanger in dollars 


SPeclite Neat at, constant pressure of heating 
sO ILL. abigy Mishel Way as 


opecific heat at constant pressure of working 
fluid in Btu/lb_ Ce 


Cost of materials for shell and tubes in dollars 


HanricavdomecOses ims cOllars 
aie Cal LONE COS Se Del cules nm doi) lars 
inner diameter of tubes imerc. 

Outer dianevet Gis cubece ih tc. 
Frictional pressure dron in lb/ft?. 
neae exchanger el tectiveness 


Mean temnverature difference ratio 


Mass flow velocity of heating fluid in 1b,/hr. 


Magis Lloweveloerty sor Wworne ei Mud in bb an. 


Heavy Transfer vecoettitcrtenry set Mmeating fluid an 
Bew/nr., ft aoe 


reat transfer coefficient of working fluid in 
SiO Mopermml Gs: |<) 5. 


Thermal “Conductivity foe Neaetine fluid in 
ny Aobe eee 5. lm 


Tiermnakveonauctiviiewm of Vertis L£iluid in 
orm bW Ag hanmaskin, 203 


Length of heat exchanrer in ft. 
MUMDeCr Of  TApUds Gatawr ous 
Mass flow rate of heating fiuid in 1b,/hr. 


Loge mean ternerature difference in °F 
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mmol OF VARIABLES ({(E.2) (cont). 


MFRYW 
MUH 
MUW 
N 
PHI 
PH2 
PY 
PH2 
PRH 


PRW 


REH 
Ri 
ROH 


ROW 


ca 


THC 
THH 
TVC 
TWH 
TLF 


TMP 


Mess (ley aiieOn \71O x: raeine iee lems <7 iar: 
NEES COs 3G) Ou sea ule i ieliened trie inte ft G 
WIM SiGle ah ie/ oj Kena gan TILOELCl DL gl a Ao Apbg as Gia. 
NUN a On mics 

[met oressUrerote heating will ine tp it. 
Gutlec pressure wot seat Mioehl Uterine ib, ata. 
ll guiyeae onwatois| Ducts) (eye (plone veal quer abdby alia) @alje) eiljoy//\ine oe 
Curlet Nressure er vOrkine fluid sins ib/ie-. 
Prandtl] number of heating wwuld 

Prandtl number of working fluid 

Neat transfer rate in Btu/hr. 

Reynolds number of heating fluid 

Reynolds number of working fluid 

Devas? one Joveereseye IIS Sia Winyeyde 
Density Pom storie dG bf Gere 

Tube svacing 

Cooler temperature of heating fluid in oF 
hotter temperature or featine fluid an CF 
Cooler temnerature of working fluid in °F 
Houter tenperature or woruing fluid in OF 
Tube Lenevihsiac tor 

BUDe) Maven ial, accor 

Overall heat transfer coefficient in Btu/hr ft? OF 


Temperature coefficient 


Sth fe 
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Appendix F. CROSS-FLOW COOLER HEAT EXCHANGER PROGRAM 


F.1 Instructions for Oneration and Inputs 

Written in FORTRAN IV, this program is designed to be run 
on the Interdata computer of the Joint Civil and Mechanical 
Engineering Computer Facility at M.I.T. With modifications in 
mae control, input and output cards, this vrogram can be run 
on any computer that is designed to comnile FORTRAN IV. 

The first three cards of the input deck contain propver- 
ties for the fluids. The next ten cards are used to input the 
mim by ten matrix F. A single number, M, is placed on the four- 
teenth card. This is the number of geometric arrangements that 
will be examined. The invut cards that follow are in two card 
sets. The first of these two cards contains an identification 
name. The second contains the data on the geometrical arrange- 
ment of the exchanger, DOUT, DIN, S, N, and KMAT. With this in- 
muy arrangement, M different geometric arrangements can be ex-— 
amined using tne Samewsriuid provertacs, 

Mie 1OUL data 1s punched on cards as follows: 

(1) All values will be in units of l1b,, hr., 
feeter Stu. and oF 

(2) All values excevt M will have fixed point 
numbers and will have a decimal point. 
Mis a floating point number 


(3) No commas are needed to senarate the 
Pape slate Onlmaesd nme card 


(4) The input values must be in the designated 


column but need not be right or left justi- 
fied 


=a 





(5) All fixed point numbers can have up to 
4 (four) digits after the decimal point. 


(6) The fluid properties are average values 
averaged over the temperature range 


CARD COLUMN 
Teele W642 arse susie Lg) Gla ye 
a: MFRW PCH Ce TWH TWC PCl 
2 PC2 PW1 PYW2 CPW CPC MUC 
3 MUW KW KC ROW ROC 


Cards 4 through 13 contain values of the matrix F 


These values are vunched 10 to a card with 6 columns per value 


1---4 
14 M 
5 -------- 18 
5 Identification name 
1--12 13--24 25--36 37--48 49--60 
G6 - DIN DOUT S N KMAT 
ill 7 Same as 15 
18 Same as 16 


These last two cards are repeated M times 


Tiesoucput Will anpear with the identiitecation mame at the 


map of a Dace and the following orinted below it: 


MASS FLOW RATE WORKING FLUID 


MASS FLOW RATE 
NUMBER OF TUBES 
INNER DIAMETER 


( 


COOLING LUI 


OF TURES 


cont. overleaf) 


=oo= 





OUinn DIATE OF SUBES 
SEACING Gf -1UbES 

EENGLA OF TUBES 

PROP 2H PRESSURE GORKIEG FLUID 
DROP fi PRESSURE, COOLENG  2huLp 
COMPARITIVE COSTS 


memeeach identification name there will be a separate nage. 
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fee uelol OF VARIABLES 


AREA Heat transfer surface area in ft?. 

COST VOst of exchnancer in dollars 

CPC Speci Pres neat sae ecOonStanuespressure of “cooling 
a) 


eM EHG! Gol Ee by ae) °F 


CPW PPeciiae Neat a@ueconsuany Dressure of working 
PiusdS in Beus lb, OF 


CTB Cost of materials for shell’ and tubes in dollars 
CW HabricationecOsts In dollars 
CWV MabEereatilon Costs oer ctube-2n dollars 
DIN Hiner sdlamever some Tubes inet. 
HOUT Cucer sdLaMnecenwsotsunoesuan ft. 
DLPC M¥LeCLionalepressure dropeor cooling fluid in 
IAs eal 
DTLM Log mean temperature difference in °F 
EC Heat exchanger effectiveness 
iF Mean temnerature difference factor 
Ge | Mass flow velocity of cooling fluid in lb,/hr. ft?. 
GW Mass flow velocity of working fluid in lb,/hr. ft?. 
HC Heau transfer eecefticien> of Cooling fluid in 


Be bv Agree ala re 222) 


HW Peart transtervcoelif etemuson working Fluid in 
Be Wine oe Tit eme ore 


KC Trertmaleconaieciv1 ty otecool~ne fluid in 
Denby en core Riv meme ce 


KV] ipernal seCOnductl Vite Of Workme Tluid in 
IS GU Me areme te em) 


LENGTH Leneth or Neauvexchancer im, it. 


ee 





ie 2 LIST OF VARIABLES (cont). 


M Number of input data groups 

MIRC Mass flow rate of cooling fluid in 1b,/hr. 
MFRW Mass flow rate of working fluid in lb,)/hr. 
MUC Viococw ty Otmeool | Mew flu Gustiae rial. 

MUW WiScos. ty Ofewomelne f1utid ne lb/ft. 

N Number of tubes 

PCl iikewepressuneson scooting filusdeine lb, fc. 
PC2 Oubletipmessiretomecooldnem mini dane iby ute. 
PW1 Inlet pressure of working fluid in lb/ft?*. 
PW2 Outlet pressure of working fluid in 1b/ft°. 
PRC Erendy rl nunbver “of ,cooling 7 luid 

BR Prandtl number of working fluid 

Q Heat transfer rate in Btu/hr. 

REC Reynolds number of cooling fluid 

REW Reynolds number of working fluid 

ROC Density of cooling: tluideimeib-/ft*. 

ROW Bensity Of Working st lid gine lb/ft ~. 

5 Tube spacing 

TCC Cooler temperature of cooling fluid in Op 
TCH Hotter temmerature of cooling fluid in °F 
TWC Cooler temverature of working fluid in OF 
TWH Hotter temperature of working fluid in OF 
LE Tube lenmen factor 
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Ps c iets On VARTABLES S)Ceomt:). 


TMF Tube material factor 
U Overall heat transfer coefficient in Btu/hr. ft?. OF 
Z Temperature coefficient 
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Appendix G - EXAMINATION OF THE ZENER CYCLE 


ies aL. Examination 

To examine the usefulness of the comouter vorogrmas de- 
veloped, the thermal cycle provosed by Prof. Clarence Zener(1) 
was examined. The vorogram for designing condensers and the 
program for evaporator design were used. 

In his pvrovosal, eee suggests the use of ammonia as 
mae working fluid. The condensor and the evaporator for this 
cycle were designed using ammonia and freon-21 as working fluids. 
Because of its high latent heat of vaporization about 540 Btu/1b,, 
ammonia requires a much larger heat exchanger for evavoration 
Pietedocs t1reon—-21 with a latent heat of vaporization of 102 Btu/7lb_. 
This difference in size can be seen in figures G-l, G-2 and G-3. 

The proverties of Zener's vorovosed cycle, as given In Ref. 1, 
Mere UuSeG aS induts to the computer programs. The allowable 
pressure adrov over the length of the condenser and evavorator 
was used as a constraint on the size of the exchangers. An 
inner diameter of % inch was used as a starting diameter. If 
this diameter caused a frictional oressure drop greater than 
the allowed vressure dron, the diameter of the tubes was in- 
creased and the heat exchanrer redesigned. 

In the condenser design, two values of the coefficient of 
heat transfer for condensation were used. The vrogram was first 


run with the coefficient being calculated as discussed in section 





a 





2.3. The design was then revised using a coefficient of 
800 Btu/hr-ft? - °F. This is the value thatGregorig(33) ob- 
tained using fluted tubes. This is further discussed in sec- 
meOnmeeso. cener in his proposal Suggested the use of this tyve 
Cie scuriace to increase the heat transfer coefficient. 

The two programs were run with the variation mentioned 
above and with different geometries.- The size of the heat ex- 
changers required for these different arrangements were then 


calculated. 


G.2 Results 

The conaensor was sized using different numbers of tubes. 
The resultant required length of heat exchanger for each differ- 
ent number of tubes are vlotted in Figure G-1l. As there is a 
significant improvement in reducing the size of the required 
condensor with the use of Gregorig surfaces, it is that size 
that is vlotted. Plotted on Figure G-l are the results for 
both ammonia and freon-21. The minimum allowable diameter of 
the tubes for each given number of tubes is also plotted. 

Figure G-2. 

The evaporator was sized in tne same manner as the con- 
denser. Zener suggests the use of surface enhancement for the 
evaporator also. However, with the high mass flow rate surface, 
enhancement annears to be of little advantage. The required 
length for the evavorator with the use of both ammonia and 


Preon—21 is plotted in Figure G=—3. 
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G.3 Conclusions of the Zener Proposal 

The use of ammonia as a working fluid in this cycle apvears 
to be totally unrealistic. The minimum size reauirement for 
the evaporator would amount to over Chiriy. bation cubic feet 
with a length of over four thousand feet. The condenser using 
meucrea Cubes would be a great deal smaller, but still is over 
SiempalTen cubic feet in total volume. 

Freon-21 is a much detter working fluid for these heat ex- 
ehanesers. The reduction in latent heat required to vaporize 
the working fluid reduces the size reauirement of the evavorator. 
Still the total volume required for the evavorator is two billion 
cubic feet. The condenser for a freon-21 system would be over 
: three hundred million cubic feet in volume. 

The sizes of the above mentioned heat exchangers would 
be for production of 100 Mega-watts of electrical vower. This 
is not an extremely large power vlant by today's standards in 
ferms of power outout. 

From a thermodynamic point of view, this cvcle is a feasible 
one. However, when the size of the required heat exchangers is 
Pakenme into account, the cycle quickly appears to be unrealistic. 
With a temperature difference of only two degrees Centigrade, in 
BINS: evaporator, the Mass flew mate is so ereat that the cycle 


is unusable. 
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